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The flavor changing neutral current processes Kl — > 7r°e + e~, S — > A s e + e~ and 
Qlh| S — > X s /j, + fi~ are studied within the minimal supersymmetric extension of the Standard 
Model. We first examine the rates for these decay modes in the MSSM with a universal 

D 

soft supersymmetry breaking sector at a Grand Unification scale. We later relax the uni- 
versality condition and investigate the FCNC transitions in a more general class of models 
with negligible flavor violation in squark mixing matrices. We find that the MSSM predic- 
tion for the kaon channel's branching fraction differs from its Standard Model value by at 
most 30% over the entire allowed parameter space. On the other hand, supersymmetric 
contributions could potentially enhance certain B — > X s £ + £~ observables by more than 
100% relative to Standard Model expectations. The impact of supersymmetry upon the 
B meson modes is strongly correlated with the MSSM value for the Wilson coefficient of 
the magnetic moment operator that mediates B — > X s ^f. 
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1. Introduction 



The Minimal Supersymmetric Standard Model (MSSM) remains one of the few well- 
motivated extensions of the Standard Model which has survived precision electroweak 
measurements [Q. The need to subject this theory to other complementary experimen- 
tal tests has consequently grown with time. Flavor Changing Neutral Current (FCNC) 
phenomenology represents one area where data are stringently confronting the MSSM 0. 
Deviations from the Standard Model may be observed in FCNC processes long before 
superpartners are detected at high energy colliders. Alternatively, failure to detect such 
departures places constraints upon weak scale supersymmetry and all other theories of 
physics beyond the Standard Model. For example, the recent CLEO observation of inclu- 
sive B — > Xg'j decay rules out charged Higgs bosons lighter than 260 GeV in Two-Higgs 
Doublet models || . This lower bound lies far beyond the reach of present direct searches. 
The B — > Xs'j measurement similarly constrains the MSSM, but its restrictive power is 
diminished by possible cancellations between different superparticle contributions |||. It 
is therefore important to study the sensitivity of other FCNC processes to supersymmetry 
and determine their limiting capabilities. 

In this article, we investigate the rare decays K L — > n°e + e~, B — > X s e + e~ and 
B — ► X s fi + fi~ within the MSSM framework. Positive signals in these channels are expected 
to be observed within the next few years provided their rates do not lie significantly 
below Standard Model predictions. The current experimental upper bound on the first 
process Br(K L — > 7r e + e _ ) e xp < 4.3 x 10 -9 is three orders of magnitude larger than the 
anticipated Standard Model branching fraction. Yet this mode is expected to be detected at 
Fermilab following completion of the Main Injector H . Present CLEO and CDF exclusive 
limits on the second and third channels, Br(S° — > X* e + e _ )cLEO < 1-6 X 10 -5 and 
Br(S° -> A"*V i >~)cdf < 2.1 x 10" 5 [g, lie within an order of magnitude of Standard 
Model predictions. Evidence for short distance b — > s£ + £~ decay may therefore soon be 
seen with the upgraded CLEO detector and in the full Run lb Tevatron data set. 

K L — > 7r°e + e~ and B — > X s £ + £~ decay have both been extensively studied in the 
Standard Model. The impact of conventional QCD upon their rates can a priori be com- 
parable to that from any new physics. Significant theoretical effort has therefore been 
directed during the past several years towards determining the precise size of strong in- 
teraction corrections to these weak transitions |]9|-[T2|. Progress has also been made in 



estimating the hadronic matrix elements that characterize the long distance aspects of 
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these processes. Chiral perturbation theory calculations indicate that the CP conserving 
component of K L — > 7r°e + e _ decay is significantly smaller than the total CP violating 
contribution, and the direct CP violating portion is believed to dominate over its indirect 
counterpart [13, 14]]. This kaon mode can thus provide an important window onto the 



nature of CP violation. 

Much less is known about the rates at which the two semileptonic FCNC reactions pro- 
ceed within the MSSM. The sensitivity of K L — > n°e + e~ to new supersymmetric physics 
has received little attention. We therefore examine the maximal variation in the MSSM 
rate for the kaon process relative to its Standard Model value in this article. Supersym- 
metric contributions to B — > X s £ + £~ decay were previously considered by Bertolini et al. 



in ref. |I(J. These authors' conclusions need to be updated in light of the B — > X s "f mea- 
surement. The important constraint placed by the CLEO observation upon the allowed 
MSSM parameter space was incorporated into the more recent B — X s £ + £~ analysis of 



ref. ||17|| , but the impact of a universal form for soft supersymmetry breaking terms was 
not studied in this latter work. As we shall see, inclusion of both the B — > X s ^ restriction 
and the universality constraint disallows sizable deviations from the Standard Model in 
the integrated B — > X s £ + £~ decay rate. Our work thus builds upon and extends previous 
supersymmetric FCNC investigations presented in the literature. 

Our paper is organized as follows. In section 2, we review the elements of the MSSM 
which are relevant to our FCNC analysis. We then discuss current restrictions upon the 
MSSM parameter space in section 3 and describe two different procedures for mapping out 
its allowed regions. We use these scanning algorithms to examine the impact of supersym- 
metry upon K L — > 7f°e + e~ and B — > X s £ + £~ decay in sections 4 and 5. Finally, we close 
with a summary of our findings in section 6. 

2. The Minimal Supersymmetric Standard Model 

The basic structure of the MSSM is well-known and has been thoroughly discussed 
in the literature ]T8|-|21J. We therefore recall just those aspects of the theory which are 
pertinent to K L — > Tv°e + e~ and B — ► X s £ + £~ decay. We first display our nomenclature 
conventions for matter superfields and their left handed fermion and scalar components 
in table I. The fields listed in the first five rows of this table carry a generation subscript 
which ranges over three family values. They are also assigned negative parities under a 
discrete Z2 symmetry in order to forbid baryon and lepton number violating interactions. 



1 This favorable hierarchy for CP conserving and violating contributions to K L — > n°e + e~ 



decay has recently been challenged in ref. 15|. 
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Table I. MSSM matter content 



The Higgs fields appearing in the last two rows transform positively under this matter 
parity symmetry. 

The chiral superfields in table I enter into the superpotential 

W = piRxRi + YVQ % U C 3 H 2 + YjjQiDjHi + Y^L^Hx (2.1) 

which governs the supersymmetry preserving interactions among matter fields. II After 
vector superfield terms are included, the supersymmetric Lagrangian schematically appears 

2 Our sign convention for contracting two SU (2) doublets is exemplified by the expansion 
H1H2 = H1H2 — H± H2 of the superpotential jjj term. 
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in component form as 



£susy = -\f^ v F^ v + X^iIPabK + {D^)\D^) + ^ilp V 
dW\* fdW\ 1 / d 2 W , , 



(2.2) 



The index G labels the color, weak isospin and hypercharge factors in the Standard Model 
gauge group, and indices A and B range over the nonabelian subgroups' adjoint represen- 
tations. All MSSM scalars are assembled into cf>, while matter fermions and gauginos are 
respectively contained within the four-component left handed ip and A fields. 

Since supersymmetry is manifestly violated in the low energy world, the MSSM La- 
grangian is supplemented with the soft supersymmetry breaking terms 

£ soft = - - [m- g g aT C~g a + m^W^CW 1 + m~ B B T CB + h.c] - m\h\hi - m\h\h 2 

- QliMlhqj - ufiMl)^ - df(Mihd° - l\{M\) l3 i 3 - l^iMl)^ (2-3) 
+ [A u i3 qiU C jh 2 + A° j q i d c j hi + A^l^hx + B\lh x h 2 + h.c] . 

In order to cut down the number of free parameters which enter into this expression 
to a manageable size, some relations among soft sector masses and couplings must be 
adopted. We shall first assume that the weak scale values of all the parameters in ( |2~3| ) 
are simply related to GUT scale progenitors. The running gaugino masses m§(/i), m~ (//) 
and m~(n) then unify at fi = Mqut just like the gauge couplings. We also equate all 
scalar mass parameters with a single mo at the GUT scale. Finally, we set the trilinear 
interaction matrices A U ' D ' E equal to AqY u,d,e at \i = Mqut where Aq denotes a common 
proportionality constant. S Imposition of this universal structure upon soft supersymmetry 
breaking terms lends predictive power to the MSSM, but the assumed simplifications are 
quite strong. We will therefore later relax some of these constraints and investigate a more 
general class of supersymmetric models. 

Renormalization group evolution of MSSM parameters down from the unification scale 
can generate a vacuum instability [pT| -p3|. As the couplings in the scalar potential run, 
the neutral Higgs fields may at some point develop nonzero vacuum expectation values 
(hi) = vi/y/2 and (h®) = V2/V2 which break electroweak symmetry. The numerical value 



We do not assume any a priori relationship between Aq and B in (|2. 



4 



v = \/v\ J rv\ = 246 GeV for their mean is fixed by the W boson mass. But the ratio 
tan/3 = Vijv\ remains a free parameter in the model. We restrict this ratio to the range 
2 < tan (3 < 55 so that Landau poles do not develop in the top or bottom Yukawa couplings 
anywhere between the weak and GUT scales. 

Electroweak symmetry breaking induces mixing among MSSM fields. In the matter 
sector, primed mass eigenstates are related to unprimed gauge eigenstate counterparts as 
follows: 



u' = T U 



S Ul u 

v! = S u -u + S U «C^ T " ~ * \S U * " V 

„ rp _ / CD L J 

d' = S DL d- S° R Cd c d' = T D 1 



S Dr d c 

S EL e- S Er C^ t ( qE l 



v' = S L V ~l -pJV qE 



(2-4) 



e 

-S Er e~ c 



The unitary S and V transformations rotate fermion and sfermion mass matrices into real 
and diagonal forms. The 3x3 quark and lepton mass matrices are simply related to the 
Yukawa couplings in the superpotential: 

Mu = ^P s U ry ut s uJ 



V2 



ME = ^P s E R y E T s E L \ 

V2 

On the other hand, the 6x6 squared mass matrices for the squarks and sleptons look much 
more complicated and involve many parameters from both the supersymmetry conserving 
and violating Lagrangians in ( |2.2| ) and (|2.3| ): 



M- = 

u 



( ^M|^ t +M2 + ^(3-4sin 2 ^)cos2/? flM v cot (3 - ^^. S Ul A u *S u ^ \ 

6 V 2 

IL*M V cot/3 - v J^ S Ur A uT S u ^ S Ur M? c t S Ur ^ + M 2 + sin 2 9 cos 2/3 

y/2 3 / 



a 



( S D ^M?S D ^ + M 2 - ^-(3- 2sin 2 9) cos2/3 LLM D tan/3 - V -^- S D ^ A D * S Dr] \ 

q ' 6 V2 

H*M D ttmf3- V -^-S DR A DT S D ^ S Dr M% t S Dr] +M 2 D - ™* sin 2 #cos 2(3 , 

V2 d 3 / 



pr>t 
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t , 1 



Mj = F N (S El M~S El + -m 2 z cos 2(3) Y 



TAf t 



Mi = 

( S E ^M 2 S E ^ +M 2 E - ^i(l-2sin 2 #)cos2/3 fJLM E tan/3 - ^5 £ M £ *S Bst 

2 ^2 

fjj* M E tan/3 - !i^^HA ET ^ t S Er M~ c t S Er ^ + M 2 E - m\ sin 2 cos'. 

2 (2-6) 
Mixing also takes place in the gaugino and Higgs sectors. The physical Dirac chargino 
and Majorana neutralino eigenstates are linear combinations of left handed Winos, Binos 
and Higgsinos: 

W~ 



X 



u 



+ V*C\ 



A M 



N 



W 3 
hi 
\h°J 



+ N*C 



(B_ t \ 

w 3 T 

— — T 

The unitary transformations U, V and iV diagonalize these fields' mass matrices 

x \ \[2m w cos (3 —\1 J 

and 

/ m B — m z sin6>cos/3 m z sin6>sin/3 \ 

AT 



(2.7) 



rriw m z cos 9 cos (3 —m z cos 9 sin (3 

m z sin 9 cos /3 m z cos 9 cos /3 fl 

\ m z sin 6* sin /3 — m z cos6>sin/3 /i 



iVt. 
(2.9) 



(2.10) 



Similarly, charged scalar mass eigenstates are combinations of hf and h^' 

f n ± \ _ f cos/3 sin/3 \ f hf 
\h ± J ~ V- sin/3 cos/3 / \hf 

The 7T 111 would-be Goldstone bosons are absorbed via the Higgs mechanism into the lon- 
gitudinal components of the gauge fields. But the remaining bosons represent 
genuine propagating scalar degrees of freedom whose tree level squared masses equal 
m 2 h± = m 2 w + m 2 +ml + 2\fl\ 2 . 

After the gauge eigenstate fields in the supersymmetric Lagrangian ( |2.2| ) are rewritten 
in terms of their mass eigenstate counterparts, @ it is straightforward to work out the 



We suppress primes on mass eigenstate fields from here on. 
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interactions of gluinos, charginos and neutralinos with quarks and squarks. We list below 
the resulting terms which participate at one-loop order in di — > dj£ + £~ decay: 

8 

C 9A = " V2g^YM^^ DljP - ~ ? DR P+)T a d 



0=1 

4 

1=1 1=1 



+ ^X 7 u\X?-P_ +Xy«P + )d + J2(xL) I dHZ? L P- +Z?*P+)d + h.c. 

(2.11) 



where 



(2.12) 



xY* = 92 [-v i \t Ul + v; 2 r u « M \ ] k 

V2m w smp 
V lm w cos p 

Z?» = - 4[| tanfliW- + iv, 3 r^ M ° - ]. 
V2 3 m w cos p 

Flavor mixing enters into these interactions through the Kobayashi-Maskawa matrix K = 

S Ul S Dl ^ and the 6x3 block components of T u and T D : 

pf _ I -pU L -nU R \ 
1 6x6 _ I 1 6x3 1 6x3 J 

2.13 

pfl _ ( yD L y D R \ 
1 6x6 ~~ I 1 6x3 1 6x3 J ■ 

Other gauge boson and Higgs terms which mediate the FCNC processes of interest are 
similarly extracted from the Lagrangian. The Feynman rules for all these interactions may 



be found in the literature 18,24 



Having set up the basic MSSM framework, we are now ready to explore its large 
parameter space. We take up this topic in the following section. 



3. MSSM parameter space 

Before predictions can be derived from the Minimal Supersymmetric Standard Model, 
explicit values for the parameters in the superpotential (|2J]) and soft supersymmetry break- 
ing Lagrangian (|2.3[ ) must be specified. In order to reduce the size of the parameter space, 
we initially adopt the common assumption that MSSM masses and couplings are simply 
related at [i = Mqut- This ansatz is motivated by the simplest supergravity theories 
25[| . A universal soft supersymmetry breaking sector at ji = Mqut also decreases the 
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likelihood of generating unacceptably large FCNC amplitudes at the weak scale. Instead, 
the magnitudes for such amplitudes in the MSSM are anticipated to be of the same order 
as those in the Standard Model. Experimentally discriminating between the two theories' 
predictions for various FCNC transitions rates will not be an easy task. But this goal is 
hoped to be achieved by a number of experimental programs within the next several years. 

In our analysis, we take as input parameters the dimensionful soft sector quantities 
Ao(Mqut), ^o(-^gut) and m^(m z ), the dimensionless ratio tan/3, and all Standard 
Model fermion and gauge boson masses and couplings. We also restrict the source of CP 
violation in the MSSM to stem from just a single phase in the KM matrix. Imaginary parts 
of and BjJj can be rotated away by field redefinitions, but Aq and \1 generally remain 
complex. The phases of these last two parameters are tightly constrained by neutron 



electric dipole moment limits [ZS]. We shall simply take Aq and fl to be real. 

In order to determine the numerical values for all the couplings in the MSSM, we follow 
a lengthy yet straightforward procedure. We first locate the GUT scale Mqut ~ 10 16 GeV 
by evolving the SU{2) L and U(1) Y gauge couplings up to the point where they meet. We 
then choose specific values for A (Mgut), ^o(-^gut), m ^( m z) and tan/3. A large value 
for the trilinear scalar coupling in conjunction with small values for the common scalar 
and gaugino mass parameters tends to yield stop masses which are too light to satisfy 
direct search constraints. Large values for either mo or m~ lead to squark decoupling 
and negligible supersymmetric contributions to FCNC decays. We therefore restrict the 
magnitudes of the three dimensionful quantities to be less than 1 TeV. 

After the gaugino, scalar and Yukawa terms in the soft supersymmetry breaking sector 
are evaluated at fi = Mqut and run down to \i = m z , the numerical values for all MSSM 
parameters except \1 and B are determined. The tree level relations 

2 mj^ sin 2 /3 — m 2 cos 2 /3 1 2 
~ cos 2/3 2 mz ( 3 ^ 

Bfl = - sin2/3(m 2 + m\ + 2\fl\ 2 ) 

fix these last two quantities up to a twofold ambiguity in sgn(Zi). Points in the MSSM 
parameter space which yield negative values for \/Jj\ 2 or B/J, fail to break electroweak 
symmetry and are rejected. Necessary tree level conditions for the existence of a stable 
scalar potential minimum 

\B/J J \ 2 >(mi + \/J J \ 2 )(m 2 2 + \IJ J \ 2 ) 

(3.2) 

(m 2 - mj) cos 2(3 > 
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must also be satisfied. We further require that the MSSM particle spectrum be consistent 
with present limits from direct superpartner searches 0. In particular, we impose the 
recent LEP 1.5 lower bound of 65 GeV on the chargino mass |[27|| . 

The final constraint which we place upon the MSSM parameter space comes from 
B — > Xg'y decay. Supersymmetry modifies the Standard Model prediction for the rare 
radiative rate by adding extra contributions to the Wilson coefficients C 7 ( tti w ^ and C%(jn w ^) 
of the electromagnetic and chromomagnetic moment operators in the AS = 1 effective 
Hamiltonian f2£|. We neglect Cs since it accounts for only 3% of the Standard Model 
b — > 57 amplitude and is not expected to be significantly more important in the MSSM. 
But the charged Higgs, chargino, neutralino and gluino contributions to C 7 which are 
tabulated in Appendix A can be quite substantial. We therefore calculate the ratio 

C 7 (m w ) M ssM fo Q x 

Kl = ~FT~i S \ 6 - 6 > 

at each point in the MSSM parameter space. We throw away all points whose values for 
R 7 do not lie within the allowed intervals 

0.4 < # 7 < 1.2 or -4.2 < R 7 < -2.4 (3.4) 



that take into account current experimental errors [[| and theoretical uncertainties [2J 

After scanning over the MSSM parameter space and imposing all the above criteria, 
we identify a number of general features which hold everywhere in the allowed regions 
except in the very large tan (3 domain: 

(i) Sizable supersymmetry contributions to K L — > 7r°e + e~, B — > X s £ + £~ and B — »■ X s j 
decay mainly arise from charged Higgs and chargino exchange. 

(ii) Flavor violating entries in F u hardly affect these FCNC processes. 

(iii) The first two generations of up and down squarks are almost degenerate. The first 
two generations of sleptons and sneutrinos are also nearly degenerate. 

(iv) Left-right squark and slepton mixing is negligible for the first two generations. 



Similar observations have previously been noted in ref. [16]. 

These characteristics provide useful guidelines for establishing less restrictive con- 
straints on supersymmetric extensions of the Standard Model than those which underlie 
the MSSM with a universal soft breaking sector. Rather than starting with a unified set 
of GUT scale parameters and evolving them down to low energies, we can instead survey 
all possible values for MSSM couplings and masses at the weak scale for which conditions 
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(i) - (iv) are satisfied. This alternate mapping procedure provides a useful check on the 
sensitivity of FCNC results upon the assumed form of the soft supersymmetry breaking 
sector. The particular weak scale quantities which must be specified in order to determine 
supersymmetric contributions to di — > dj£ + £~ decay are listed below: 

• a common mass for the superpartners of left handed up and charm quarks, 

• the masses mi L , ttii r and mixing angle ai for top squarks, 

• a common mass for the first two generation sneutrinos, 

• the Wino and charged Higgs masses and m^± , 

• the superpotential \1 parameter and tan/3. 



This parameterization is similar in spirit to the one adopted in ref. |T7j . As in our universal 
soft sector analysis, we shall restrict the dimensionful quantities to the sub-TeV regime 
and restrict the dimensionless VEV ratio to 2 < tan (3 < 55. 

In the next two sections, we will investigate the K L — > 7r°e + e _ and B — > X s £ + £~ 
transitions utilizing the universal soft sector scanning procedure as well as the more general 
mapping method. The greater labor required to explore the much larger nine-dimensional 
parameter space in the latter approach is offset by several simplifications. For example, 
neutralino and gluino contributions to di — > dj£ + £~ may simply be neglected. The T u 
matrix also reduces to the nearly diagonal form 



(I 
1 




Vo o 






cos ax 



sin ax 









1 
1 







— sin a- t 



COS Off / 



(3.5) 



when criteria (i) - (iv) are imposed. But most importantly, no renormalization group 
evolution needs to be performed. Searching the nine-dimensional parameter space for sets 
of points where supersymmetry significantly enhances or suppresses the rare decay modes 
is thus rendered tractable. 



4. K L — ► n°e + e decay 

The total amplitude for K L — > 7r°e + e _ decay can be decomposed into CP conserving 
and violating parts. The former starts at second order in the electromagnetic interaction. 
As a result, the CP conserving branching fraction Br(K L — > 7r°e + e~)cp — (0.3-1.8) x 10 -12 
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is significantly smaller than its CP violating counterpart | H%[ . Moreover, the indi- 
rect component of the CP violating amplitude is believed to be smaller than the di- 
rect part. Present data imply Br(K L — > 7r°e + e~ indirect < 1-6 x 10~ 12 |T|], while the 



Standard Model prediction for the direct CP violating contribution lies in the range 
Br(i^ 7r°e + e-)direct ^ (2.5-9.0) x 10" 12 JTIJ. § The CP conserving and indirect CP 
violating amplitudes may be computed using experimentally determined values for chiral 
Lagrangian coefficients and the e parameter which automatically include possible MSSM 
contributions. Since these quantities are consistent with Standard Model predictions, we 
will focus exclusively upon the direct CP violating component in our supersymmetric 
FCNC investigation. 

The analysis of K L —>■ Tr°e + e~ decay is greatly facilitated by working within an effec- 
tive field theory framework. In this approach, heavy degrees of freedom are successively 
integrated out from a specified full theory, and the resulting effective theory is run down 
to a low energy hadronic scale using the renormalization group. A finite number of flavor 
changing dimension-6 operators generated by this process at /i = m w enter into the effec- 
tive Hamiltonian 7i er r which governs the dynamics of the low energy theory. The complete 
set of left handed AS = 1 four-fermion operators that originate from the Standard Model 
and mediate K L — > 7r°e + e _ are catalogued in ref. [II |. If the starting full theory is taken 
to be the MSSM, additional terms with right handed flavor changing currents appear in 
Heft- But since their coefficients are tiny, the extra operators may be neglected without 
loss. 

After the effective theory is evolved to low energies, the direct CP violating 
K L — > 7r°e + e~ amplitude is well approximated by the long distance matrix element 
{7T°e + e~ \H e a\Ki) of the truncated Hamiltonian 

—G 

H eS = —^-K* ts K td [yrvQrv + ViaQia] + h.c. (4.1) 

All information associated with short distance physics is encoded into the Wilson coeffi- 
cients 



V7V 



2tt sin 2 9 

OLem 

Via - 



2n sin 2 6 



Po + (Y(xt) + F SUSY ) - 4 sin 2 9 (Z(x t ) + Z SUSY ) + P E (E(x t ) + £ SUSY ) 

Y( Xt ) + y SUSY 



(4.2) 



5 The sizable uncertainty in the direct CP violating branching fraction primarily stems from 
KM angles which are poorly constrained at present. The predictions' precision should substantially 
improve when the KM unitarity triangle is better determined in upcoming B-factory studies. 
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of the semileptonic operators 

Q 7V = s^(l--f 5 )de-f fI e 



(4.3) 

Qia = 57^(1 -7 5 )de 7/i 7 5 e. 



The functions 



Y( Xt ) = 4Xt X \ + . 3X ' , 9 logs* 
1 t; 8(1 -i t ) 8 1 -x t 2 8 



108x t - 259x? + 163x? - 18xf -8 + 50x t - 63x 2 - 6xf + 24xf , 
Z( *'> = 144(1 + 72(1 -V ° gX< ' ' 



„, . 18x t -llx 2 -x? 4-16x4 + 90;? 

12(1 -x t y 6(1 -x t ) 4 

of the variable xt = (mt/m w ) 2 summarize the high energy contributions to K L — > n°e + e~ 
decay that are common to both the Standard Model and MSSM. Next-to-leading order 
QCD corrections which are the same in both theories are incorporated into Pq ~ 0.7 and 
Pe — —0.01 Therefore, all effects of supersymmetry upon the rare weak transition 



reside within the y SUSY ? j^susy an( j ^susy p arame ters in eqn. 

The values of the yrv and in a coefficients are obtained after matching renormalized 
s — > de + e~ amplitudes calculated in both the full and effective theories. We perform this 
matching at \i = m w in the Standard Model as well as in its minimal supersymmetric 
extension. The one-loop bubble, penguin and box graphs which must be evaluated on the 
MSSM side of the matching condition are displayed in figs, la, lb and lc. After a long but 
straightforward computation, we find the charged Higgs, chargino, neutralino and gluino 
contributions listed in Appendix B to the y SUSY and Z SUSY terms in ( |4.2| ). 

Several points about our supersymmetric matching results should be noted. Firstly, 
we have not included any bubble or penguin graphs which involve neutral Higgs boson 
exchange in fig. la or fig. lb. Such diagrams are proportional to the electron mass m e and 
are negligibly small. Similarly, the first box graph in fig. lc with a charged Higgs running 
around the loop vanishes in the m e — > limit and may be safely ignored. Secondly, we have 
intentionally not calculated _E SUSY which should be comparable in size to E(xt) — 0.25. 
Since the numerical value for Pe is almost two orders of magnitude smaller than that for Pq, 
the term proportional to Pe in ( |4.2| ) is negligible and _E SUSY is unimportant. Thirdly, we 
have ignored gluino loop matching contributions to charged current four-quark operators in 
the AS = 1 effective Hamiltonian. Such contributions affect yrv an d Via at the sub-percent 
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level. Finally, the matching condition results displayed in Appendix B are independent of 
any choice for the structure of soft supersymmetry breaking terms. 

The short distance y-jy and Via Wilson coefficients enter into the direct CP violating 
K L — ► 7r°e + e _ partial width along with a long distance hadronic matrix element. After 
neglecting tiny imaginary components in the former quantities and relating the latter to 
the measured K + — ► 7r e + u rate via an isospin rotation, we find [III] 



Br(K L -> 7r°e+e-) direct = Br(K+ -> 7r°eV 



r(K L ) 
r(K+) 



K,. 



2 



(yrv + yh)- (4-5) 



The sensitive dependence of this expression upon poorly known KM angles can be removed 
by normalizing the MSSM branching fraction to its Standard Model analogue. The ratio 

p(K L - *°e+e-) = MKl - *° e l e 2 MSSM = B : {Kl - *° e+e ~^ (4.6) 

Bi{K L -> 7T u e+e-) S M 

thus cleanly quantifies supersymmetric enhancement or suppression of the rare kaon mode's 
rate. 

Following the universal soft sector scanning procedure outlined in section 3, we eval- 
uate p(K L — > 7r°e + e~) in the allowed regions of MSSM parameter space. Representative 
results for two-dimensional slices through this space are illustrated in fig. 2a and fig. 2b. In 
the first LEGO block figure, p(K L — > Tv°e + e~) is plotted as a function of Ao(Mgut) and 
mo(McuT) with m^(m z ) = 90 GeV and tan (3 = 5.0 held fixed. In the second figure, we 
display p(K L — ► 7r°e + e~) as a function of m^{m z ) and tan/3 with A (Mgut) = —500 GeV 
and mo(McuT) = 100 GeV. We have taken the sign of the fl parameter in the superpoten- 
tial to be positive in both plots. Points in the MSSM parameter space that are excluded 
by one or more of our imposed criteria are indicated in these figures by LEGO blocks 
which are saturated at their maximum values. Looking at the results in figs. 2a and 2b, 
we see that supersymmetric effects in the MSSM with a universal soft sector reduce the 
direct CP violating K L — > 7r°e + e~ branching fraction relative to the expected Standard 
Model rate by at most 10%. These discrepancies between the Standard Model and its 
minimal supersymmetric extension are unfortunately too small to be realistically detected 
by upcoming experiments within the next several years. 

It is interesting to examine whether larger deviations could result if some of the 
stringent assumptions underlying the universal soft sector MSSM are relaxed. We have 
therefore performed several scans involving hundreds of thousands of points over the more 
general nine-dimensional parameter space discussed in section 3. These scans reveal that 
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pockets in the larger space exist where discrepancies between the Standard Model and 
the MSSM are three times larger in both the positive and negative directions than those 
we previously uncovered in our more restricted searches. But even a 30% supersymmetric 
enhancement or suppression of the K L — ► ir°e + e~ rate relative to its Standard Model value 
is unlikely to be experimentally resolvable in the near future. 

The potential impact of supersymmetry upon the rare kaon mode is thus disappoint- 
ingly small. But as we shall see in the next section, the prospects for detecting signs of 
supersymmetry in B — > X s £ + £~ decay are brighter. 



5. B -> X,e + e~ and B 



Inclusive B 



X. q e 



X a fji + fM decay 

and B — > X s [i + [j,~ decay share several similarities with the 



K L — > 7T e + e~ transition. Like their kaon analogue, these B meson reactions are most 

conveniently analyzed within a low energy effective theory framework. The complete list 

of dimension-6 operators in the AB = 1 effective Hamiltonian that participate in B — > 

X s £ + £~ decay may be found in ref. J29|. The Wilson coefficients of the bottom sector 

operators evaluated at the W scale are trivially related to their strange sector analogues. 

The next-to-leading order evolution of these coefficients from \x = m w to fx = m& has been 

calculated in ref. []T^[] . The details of this strong interaction running are quite complicated, 

and we will not present them here. Instead, we will simply apply the main results in our 

study of supersymmetric effects upon the rare decay modes. 

The inclusive rate for the meson level process B — > X s £ + £~ may be approximated 

by the rate for the free quark transition b — > s£ + £~ [|30[]. Two independent variables are 

required to describe the latter process. We choose them to be the rescaled lepton energies 

2Eg+ 2E^- 

y + = and y- = (5.1) 

m b m b 

measured in the b quark rest frame. When expressed in terms of these variables along with 
the rescaled squared invariant mass s = y+ + y- — 1 of the lepton pair, the differential 
decay rate looks like 



d 2 V(b^s£+£-) G 2 F ml\ K* ts K tb \ 2 / a EM \ 2 



dy+ dy. 
4 

+ - 



8(1 



16n 3 V 4tc 
S) + (l-y + ) 2 + (l 



y- 



y+i 1 - y+) + y-i 1 - y- 

(Cf) 2 



\a 



eff, 



+ Cio) 



2 + 2mj 



sm 2 J 



+ 4(1 - s)Cf Re(C 9 ett (s)) + 2(y+ - y_)C 10 2Cf + sRe(C 9 ett (s)) 



-»eff 



■*eflr 



(5.2) 
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The quantities 



,16/23 



(5.3a) 



i=i 



/~<eff 
°9 



7T 



+^) (-0.1875 +x;^ ai+i ) 



a s (m w ) rj 

+ ( 0+ 2fl < Z &) + ZSVSY ) + Wxt) + £ SUSY ) (0.1405 + V q^ +1 ) 

sin p z — ' 



Cio 



+ 1.2468 + rf- \n + Sit] + £;/i(— , s) + uMl, s) + vMO, s) 
^— ' L nih 

Y(x t )+Y i 



i=i 

-SUSY 



sin 2 # 



(5.36) 



(5.3c) 



that enter into the partial width depend upon the strong interaction coupling ratio 
rj = a s (m w )/a s (mb) and various matching condition functions which we previously en- 
countered in our K L — > n°e + e~ analysis. The effective coefficients also involve the com- 
ponents of the following 8-dimensional vectors: 

(0.6087, 0.6957, 0.2609, -0.5217, 0.4086, -0.4230, -0.8994, 0.1456), 

(2.2996, -1.0880, -0.4286, -0.0714, -0.6494, -0.0380, -0.0186, -0.0057), 

(0, 0, -0.3941, 0.2424, 0.0433, 0.1384, 0.1648, -0.0073), 

(0, 0, 0, 0, 0.0318, 0.0918, -0.2700, 0.0059), 

(0, 0, 0.8331, -0.1219, -0.1642, 0.0793, -0.0451, -0.1638), (5.4) 
(0, 0, -0.2009, -0.3579, 0.0490, -0.3616, -0.3554, 0.0072), 
(0, 0, 1.7143, -0.6667, 0.1658, -0.2407, -0.0717, 0.0990), 
(0, 0, 0.2857, 0, -0.2559, 0.0083, 0.0180, -0.0562), 
(0, 0, 0.1429, 0.1667, -0.1731, -0.1120, -0.0178, -0.0067). 
Finally, the functions h(z, s) and u(s) which appear in eqn. ( |5.36p are given by 

for x = 4z 2 / s < 1 



m 

hi 
Pi 
Qi 

Ti 
Si 
U 
Ui 
Vi 



h(z, s) 



x\ 



log 



2 2 
r Li 2 (s) - - log(s) log(l - a) - -7T 2 - 
1 3 9 3(1 + 2s 

2s(l + s)(l-2s), 5 + 9s -6s 2 

- 3(1-^(1 + 21) l0g(s) + — 



2 arctan(l/i/ 
5 + 4s 



27T, 



X 

log(l - s 



1), for x = 4z 2 /s > 1 



6(l-s)(l + 2s) 



(5.5) 
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In order to consistently compute the differential rate to one-loop order accuracy, we only 
retain terms in d 2 T / 1 dy +dy _ up to linear order in lu(s). We also set w(s) to zero in the 
interference terms in eqn. ( |5.2| ) which are proportional to Re(C| ff (s)). 

The Standard Model prediction for the B — > X s £ + £~ decay rate is simply recovered 
from these formulae by setting y SUSY ? ^susy an( j ^susy to 

zero and equating Ci(m w ^ 

and Cs(m w ) with their Standard Model values. In this case, the partially integrated rate 
dT/ds differs from the corresponding result discussed in ref. Jl^] by a term proportional 
to the lepton mass. This term was previously neglected because its effect is very small 
over nearly all of phase space. However, its contribution to the integrated decay rate is 
not suppressed by mi provided no lower cut on s is imposed. As no such cut has been 
performed in the recent CDF analysis of B — > K*/j, + fi~ ||, we retain this additional term 
in dT/ds. 

The matching condition expressions for y SUSY and Z SUSY that enter into the differ- 
ential rate for B — > X s e + e~ trivially differ from those for K L — > n°e + e~ by just the flavor 
label renamings specified in Appendix B. As in our previous K L — > Tt°e + e~ matching com- 
putation, we neglect £' SUSY and set Cg(m w ) to zero. For the B — > X s fi + fi~ mode, slepton 
and sneutrino indices must also be transformed from the first to second generation. The 
numerical values for individual matching contributions to y SUSY and Z SUSY turn out to 
be almost identical for all three FCNC channels which we consider in this paper. The fla- 
vor independence of y SUSY and ^ SUSY implies that regions of the B — > X s £ + £~ spectrum 
dominated by Cg ff and C\q are no more sensitive to supersymmetry than K L — > 7i°e + e~ 
decay. Sizable discrepancies between the Standard Model and its minimal supersymmetric 
extension can therefore only arise in B —> X s £ + £~ observables which depend to a large 
extent upon Cy ff . 

We should note an interesting point regarding the supersymmetric limit of these 
matching results. The approximate cancellation between different superpartner contri- 
butions to B — > Xg'y decay has been interpreted as a manifestation of T(b — > 57) = in 
the supersymmetric limit [Q. Since the magnetic moment operator which mediates the 
radiative transition belongs to a linear multiplet, it cannot arise in a fully supersymmetric 
effective Hamiltonian [31 1. No analogous argument can be made for b — > s£ + £~ decay. 
Penguin and box diagrams generate effective four-fermion operators that form the highest 
components of vector superfields. Such .D-terms survive in the limit of exact supersymme- 
try. 
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As low statistics will most likely hinder experimental determination of the full differ- 
ential spectrum, we need to consider various integrated observables. Following the CDF 
analysis presented in ref. ||, we first integrate the rate over the lepton pair mass regions 

m £+e - E (2m £ ,2.9 GeV) U (3.3 GeV,3.6 GeV) U (3.8 GeV,4.6 GeV). (5.6) 

These disjoint intervals exclude mi+i- values for which the inclusive B — > X s £ + £~ rate 
is dominated by intermediate J/ip and ip' states. By restricting our analysis to just this 
nonresonant region, we ensure the validity of the free quark approximation to inclusive B 
meson decay. 

It is customary to reduce the uncertainties in the b — > s£ + £~ partial width by normal- 
izing it to the semileptonic rate 



r(6 - ce+u) = (l - 2as{mb) L 2 - ^)(1 - ^) 2 + 51 ) 

v ; 192tt 3 yv m 6 ; \ 3tt [ V 4 a m b J 2\ J 



(5.7) 



which is related to the measured branching ratio Br(£? — > X c e + u) = 0.104 ±0.004 [Q. The 
function 

g(z) = 1 - 8z 2 + 8z 6 - z 8 - 24z 4 log z (5.8) 

appearing in this rate formula represents a phase space suppression factor. The sensitive 
dependence of T(b — > s£ + £~) and F(b — > ce + u) upon KM angles cancels in their ratio. 
However, errors in the numerical evaluation of the B — > X s £ + £~ partial width can be 
reduced to only the 10-20% range due to uncertainties in quark masses and interference 
effects from excited charmonium states |32| . Therefore, signals of new physics beyond the 
Standard Model will be detectable only if they significantly exceed this level. 

After integrating the differential rate in eqn. ( |5.2j ), expressing the result in terms of 
the ratios 

p C 7 (m w ) MSS M „ Y(x t )+Y SVSY Z(x t ) + Z SUSY 

C 7 {m w ) S M Y(x t ) Z(x t ) 

(5.9) 

and adopting the parameter values m c = 1.3 GeV, m& = 4.7 GeV, mt = 176 GeV, 
oi s {Tn z ) = 0.118 and \Kt s K t b/ K c t,\ 2 = 0.95, we find the following nonresonant branching 
fractions: 

Bi(B -> X s e + e~) NR = 3.0 x 10" 7 [5.5 + 2.3^ + 17.6R 2 , + 3.7R 2 Z - 2.1R 7 R Y 

(5.10) 

+ lAR 7 Rz - ll.hRyRz + 4.6i? 7 + 8-l.Ry - 5.3R Z ] 
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Br(S -> A" s //+/Onr = 3.0 x 10" 7 [2.9 + 0.8i$ + 17.5i& + 3.7fl| - 2.1R 7 R Y 

(5.11) 

+ 1.4R 7 i2 z - HARyRz + 0.7R 7 + S.1R Y - 5.3R Z ]. 

The Standard Model values 7.3 x 10" 6 and 4.9 x 10" 6 for the B -> X s e+e~ and S -> 
X s [i + fi~ branching ratios are recovered by setting R 7 = Ry = Rz = 1 in these formulae. 

The second b — ► s£ + £~ observable we consider is the lepton-antilepton energy asym- 
metry 

_ N(E £ - >E £+ )- N(E e + > E £ -) 

N(E e - > E e +) + N(E i+ > E e -)' [ ' ' 

Here N(E £ - > E £ +) denotes the number of lepton pairs whose negatively charged member 
is more energetic in the B meson rest frame than its positive partner. Since A is odd under 
charge conjugation whereas Br (5 — > X s £ + £~) is even, the information about the differ- 
ential spectrum encoded into the former observable does not overlap with that contained 
within the latter. The value for A is most simply determined for charged S ± mesons which 
do not suffer from complications associated with B-B mixing. Counting only those lepton 
pairs whose invariant mass lies within the intervals specified in eqn. (|5.6Q , we find 



3 x 10 _T 

i NR = p fT3 ' v ff+ , , [1-2 - 1.0R 7 + 4.0R Y - 2.6R z ]Ry. (5.13) 
tir{B — > A s l^l Jnr 

This expression yields 7% and 10% for B — ► X s e + e~ and B — ^ X s/ u + /i~ in the Standard 
Model, respectively. 

Since deviations of Ry and Rz from unity over the allowed MSSM parameter space 
are small, supersymmetric effects in the B — > X s £ + £~ channel critically depend upon the 
ratio R 7 . In order to isolate the FCNC mode's sensitivity to this quantity, it is instructive 
to first artificially set R 7 = 1 everywhere throughout the region of MSSM parameter 
space allowed by all constraints. We then find that supersymmetric effects upon B — ► 
X s £ + £~ are quite similar to those for K L — > ir°e + e~ . Nonresonant B — > X s £ + £~ branching 
ratios are suppressed by at most 10% relative to their Standard Model values when GUT- 
scale universality is imposed, and changes in lepton asymmetries are even smaller. If the 
universal soft sector assumptions are relaxed, B meson rates and asymmetries then deviate 
from Standard Model predictions by at most 30%. 

Of course, the value for the Wilson coefficient of the magnetic moment operator does 
not coincide in most regions of MSSM parameter space with its Standard Model coun- 
terpart. When GUT-scale universality is assumed and all but B — ► X s ^ constraints are 
imposed, we find that R 7 smoothly varies between —2.5 and 3 for large values of tan/3 
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and between 1.0 and 1.3 for tan/3 ~ 2.5. This tan/3 dependence primarily stems from the 



1 / cos (3 = \/l + tan 2 (3 enhancement of the chargino interaction matrix X^ R in (|2.12|) . If 
we instead search over the nine-dimensional parameter space discussed in section 3, we find 
that R-j ranges at the factor of two level from — tan/3 to tan/3. This potentially large R-j 
variation underscores the stringent nature of the limits in ( |3.4j ) set by the CLEO B — > X s ^ 
observation. 

It is useful to separately consider two different scenarios for B — > X s £ + £~ decay which 
depend upon the sign of Rj\ 



Case 1. We accept all points in the MSSM parameter space for which 0.4 < Rj < 1.2. 
The extremal MSSM values we then find for our B — > X s £ + £~ observables are displayed 
in table II as fractions of their Standard Model counterparts. Looking at the entries in 
the table, we see that sensitivity of the nonresonant branching ratios to supersymmetry is 
fairly minimal. The reason for the slight variation can be traced to the coefficients of the 
Ry and Rz terms in ( |5.10[ ) and ( |5.11| ) which are considerably larger than those for the 
i?7 terms. Since Ry and Rz do not vary significantly from unity and the magnitude of 
i?7 is constrained by experiment, supersymmetric effects in the MSSM with a universal as 
well as more general soft sector are never highly pronounced in the nonresonant branching 
ratios. 

On the other hand, the relatively larger coefficient of the Ry term in eqn. ( |5.13| ) 
induces greater sensitivity in asymmetry observables to supersymmetric deviations from 
the Standard Model. In fig. |3], we plot the ratio of the electron asymmetry observable 
A(B — > X s e + e~)NR evaluated in the universal soft sector MSSM relative to its Stan- 
dard Model value. This ratio is displayed in the figure as a function of Ao(Mgut) 
and mo(McuT) on a two-dimensional slice through the MSSM parameter space with 
m^{m z ) = 150 GeV, tan/3 = 20 and sgn(/i) = +1 held fixed. Looking at the LEGO plot, 
we see a substantial volume in parameter space exists where supersymmetric deviations 
from the Standard Model are sizable. 
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Observable 




MSSM with universal 
soft sector 


MSSM with relaxed 
soft sector assumptions 








minimal 


lllCLjVllllGLl 


minimal 


TH £1 VI TT1 J} 1 
111 Cl^V 1 111CL1 


Br(B - 


-> X s e+e- 


")nr 


77% 


107% 


72% 


119% 


AiB- 


-> X s e+e- 


)nr 


80% 


170% 


81% 


196% 


Br(S - 


-* X s n + ij- 


~)nr 


88% 


102% 


71% 


121% 


A(B- 




")nr 


85% 


148% 


86% 


164% 



Table II. Extremal MSSM values of B -> X s £ + £~ observables 
as fractions of their Standard Model counterparts. 

Case 2. We accept all points in the MSSM parameter space for which —4.2 < Rj < —2.4. 
Scans over the parameter space of the MSSM with GUT scale universality then reveal 
that i?7 never dips below —2.5. This limiting value coincides with the lower end of the 
experimentally permissible B —>■ X s ^ partial width range. It is important to note that 
T(B X s e + e~) is strongly correlated with V(B — > X s ^) since a large short distance 
contribution to B — > X s e + e~ comes from the small dielectron mass region where the 
intermediate photon is only slightly off-shell. Maximal suppression of the radiative rate 
therefore leads to a 10% suppression of the semielectronic rate. On the other hand, F(B — > 
X s iJL + ii-) and T(B -> X s j) are anticorrelated for R 7 < due to the CfRe(Cf (a)) 
term in (|5.2|) . The semimuonic branching fraction increases by 20% when R? = —2.5. 
More importantly, sizable signals of supersymmetry can be detected in the asymmetry 
observable. *4nr increases relative to its Standard Model value by factors of 3.6 and 2.6 
for B — > X s e + e~ and B — > X s fi + fi~ when R7 ~ —2.5. 

If we again relax the universal soft sector assumptions, we find many points in the 
more general nine-dimensional parameter space where R-j drops down to the lower end 
of its experimentally allowed range. The nonresonant B — > X s e + e~ and B — > X s fi + fi~ 
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branching ratios are then enhanced by up to 90% and 110%. Lepton asymmetries are also 
enhanced by approximately a factor of 3 compared to Standard Model expectations when 
i?7 belongs to the allowed negative range, i 

The two different scenarios we have investigated for the impact of supersymmetry upon 
B — > X s £ + £~ decay clearly have different phenomenological implications. The likelihood 
that this FCNC process will display interesting evidence for supersymmetry or else usefully 
constrain the MSSM parameter space strongly depends upon the value for R7. After more 
precise experimental measurements and next-to-leading order theoretical calculations are 
completed in the near future, the allowed range for i?7 should shrink by about a factor 
of 3. It is possible that the Standard Model prediction Rj = 1 will then no longer be 
consistent with CLEO data. Such a result would represent an intriguing finding. 

At the present time, the CLEO measurement provides no means for determining 
the sign of the b —>■ s'j amplitude. However, the current 3.5a discrepancy between the 
Standard Model and LEP data for T(Z — > bo) may suggest that R-j is negative in the 
MSSM framework. Large positive MSSM corrections to T(Z — > bb) would ameliorate the 
conflict between theory and experiment. As such positive corrections are correlated with 
negative values for i?7, our second B — > X s £ + £~ scenario may be favored. In this case, 
sizable deviations in B — > X s £ + £~ branching ratios and lepton asymmetries from Standard 
Model expectations should hopefully be detected in the next few years. 

6. Conclusions 

In this article, we have studied the impact of supersymmetry upon the FCNC pro- 
cesses K L — > ir°e + e~ , B — > X s e + e~ and B — > X s (i + (i~ . We found that the rate for the 
kaon mode does not vary from Standard Model predictions by much more than 10% in 
the minimal supersymmetric extension with a universal soft breaking sector. Qualitatively 
similar results hold for the B meson nonresonant branching fraction. Since the premise 
underlying the MSSM with GUT scale universality is not necessarily realized in nature, we 
have also considered these transitions in a more general class of models in which certain 
soft sector assumptions were relaxed. We then uncovered regions in a nine-dimensional pa- 
rameter space where Br(S — > X s £ + £~)nr is significantly enhanced relative to its Standard 

6 It is interesting to note that deviations of R7 from unity in Case 2 can induce larger variations 
in r(B X,£T£~) MSSM /T(B -> X s £ + £-) S M than in F{B -» X s7 ) MSSM/T(i? — > Xj7)sM- 
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Model value. Charge conjugation odd lepton asymmetries can exhibit even larger devia- 
tions from Standard Model expectations. So signals of supersymmetry could be detected 
in B — > X s £ + £~ decay in the next few years. 

The general approach which we have followed in this paper to investigate supersym- 
metric contributions to a particular class of rare modes can be applied to several other 
interesting processes which might reveal larger discrepancies with the Standard Model. 
For example, our analysis of semileptonic d\ — > dj£ + £~ decay can readily be extended to 
its neutrino analogue di — > djvv. Similar methods can also be used to look for weak scale 
supersymmetry in K — > fj, + n~, B — > r + r~ and B-B mixing. Theoretical and experi- 
mental study of all these processes will help to constrain whatever physics lies beyond the 
Standard Model which is still waiting to be discovered. 
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Appendix A. MSSM b — > .$7 matching conditions 

We list below the VF-scale matching contributions to the coefficient C 7 of the magnetic 
moment operator in the AB = 1 effective Hamiltonian which arise from one-loop MSSM 
diagrams: 

Standard Model graphs: 



SC 7 =^h(x t ) (A.l) 



Graphs with charged Higgs loops: 



^ = \\\^ COt2 ^ A (^T-) + M"^-) } ( A - 2 ) 

6 [2mf ± m h± m h± J 



Graphs with chargino loops: 

1 6 2 2 

sr - V" V" m ™ 

6Cr 3to^^mL 



x {-(^)L(^Wi(^) + (^)L(^)x,-^/ 2 (^)} 

L 2 Vm.+ / mi Vm_ + / J 



Graphs with neutralino loops: 

1 6 4 2 

e^y 1 \ ^ \ ^ m w 



X 



Graphs with gluino loops: 



x{-(r^)L(r^) A3 /3(^) + 2(r^)L(r^) A 3^/4(^)} 

l V m~ / nib V m~ / J 



-m 2 T N /m 2 T 



(A.3) 



(A.4) 



(A.5) 



The one-loop integral functions which enter into these matching conditions are given by 

, , , -7 + 5x + 8x 2 2x- 3x 2 . 
h{x)= 6(1 -s)3 ~0^xY gX 

2(1 -x) 2 (1-x) 3 (A.6) 

/3(X)= 6(1 -xf + Jl^W l0gX 
1 + x X 

h{x) = w^w + ^w gx - 
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Appendix B. MSSM s — ► de + e and 6 — ► ae + e matching conditions 

We tabulate below the VF-scale matching contributions to the y SUSY and Z SVSY 
parameters which appear in the Wilson coefficients y-jy and y-jA of the strange sector 
semileptonic operators Q-jv and Qia in eqn. (f4.2|). These same formulae hold for Wilson 
coefficients of analogous operators in the bottom sector. The KM matrix label q and 
numerical index i respectively equal q = d, % = 1 and q = b, i = 3 for s — > de + e~ and 
b — > se + e~ decay. 

Z-penguin and bubble graphs with charged Higgs loops: 

5y SUSY =SZ SUSY = _I cot ^ Xt/5 (^L) (B.l) 

° m h± 

7-penguin and bubble graphs with charged Higgs loops: 

SY SUSY = 

1 m 2 (B.2) 

^ SUSY = -^cot 2 /3 /6 (^) 1 } 

72 m h± 

Z-penguin and bubble graphs with chargino loops: 

6 2 



^92^ts^tq A B = ll }=1 
^ Aj Aj Aj 



1 / 2 



U^ 2 >±)^ s ^/i^;i) 

^7 AJ J 



(B.3) 



7-penguin and bubble graphs with chargino loops: 
^ySUSY = 



^^sa^tt^^^^'-Mj) (B ' 4) 

Z-penguin and bubble graphs with neutralino loops: 

6 4 



D L) 

) Bi 



x {c 2 (m|o, m| A , m| B )(r^r^ f ) AB <^ - c 2 (mf v m|o, m|o )6 ab (N* 3 Nj3 - JV; 4 JV j4 ) 
- im^om^oco(m^,m|o,m|o)5 AB (A^ /3 A^; 3 - iV j4 iV* 4 )J 

(B.5) 
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7-penguin and bubble graphs with neutralino loops: 

SY SVSY = 



1 6 4 m 2 ™~o (B.6) 

r 7 SUSY_ 1 V" V" m w ( 7 Dl\\ (7,D h \ f ( *i \ K ' 



Z-penguin and bubble graphs with gluino loops: 



(B.7) 



sy svsy = ^susy = 4, ^ (r^)L(r^) BiC2 (m|,m| v m| B )(r^r^^ 

6 9 2 JX ts JX tq A B = 1 

7-penguin and bubble graphs with gluino loops: 

SY SUSY = 

„2 JL ffl 2 m ? (B.8) 

r Z SUSY_ 93 \^ m w ( V Dl\\ (r D L} f ( " L g \ ^ 1 

HI ts iq A=1 dA dA 

Chargino box graph: @ 

^ SUSY = -t^jt- E E (^)L(^U(4,^,<,<)TO , 

9%KtsKtq A ~ Xl ~^ 1 x > ^ (B.9) 

5Z SUSY = 

Neutralino box graphs: 



<^SY = 2sin 2 g ^SUSY + 2 < £ £ (^)L(^)- 

Z 92*ts^tq A=lt>J=1 

x [d 2 (m|o , m|o , , m? i ) (iV; 2 + tan 07^ ) (N j2 + tan #iV 7l ) 
+ ^m^om^odol^lo, m|o , m|^, m\ l ){N l2 + tan6W 7 i)(iV* 2 + tan6W*^ 

^ SUSY = -T^JT E £ (^)L(3?*)a< sec 2 61 

x [d 2 (m|o, m|o , m| A , mjjN^ N n + -m^m^o do(m|o, m|o , m|^, mtjiVjiiV*!] 

(B.10) 



Our chargino and neutralino box graph matching condition results differ from those reported 



in ref. [16] by overall signs. 
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The one-loop integral functions which appear within these MSSM matching conditions 
are given by 



fe(x) 
fr(x) 
fs(x) 

co(ml,7n%,ml) = 

l 2 2 2\ 

c 2 {m 1 ,m 2 ,m 3 ) 

7 / 2 2 2 2\ 

do{m 1 , m 2 , m 3 , m 4 J 



x 



+ 



X 



log a; 
+ 



l-x ' (1-x) 2 

38x - 79x 2 + 47x 3 4x - 6x 2 + 3x 4 



6(1 -x) 3 (l-x) 4 
52 - lOlx + 43x 2 6 - 9x + 2x 3 



logx 



6(1 -x) 3 
2 - 7x + llx 2 



+ 



(1 - x y 



+ 



6x 3 



(l-x)' 



(l-x) 4 
logx 



logx 



mj 



mf log 



(mf — m 2 )(mf — m 2 ' 



3 

8 



l r 

4 



2i m l 

™I log — 
Li 2 



(mf — m 2 )(mf — m 2 ) 



+ (mi <-> m 2 ) + (mi <-> m 3 )j 



+ (mi <-> m 2 ) + (mi <-> m 3 )j 



2! m l 

mf log 

(mf — m|)(mf — m|)(mf — m 2 ) 



+ (mi <-> 777.2) + (mi JTI3) + (mi <-> 7714) 



d 2 (mf,m|,m§,m|) 



4 1 TO 1 

™i log — 



(mf — m|)(mf — m 2 )(mf — m 2 ) 



+ (mi <-> m 2 ) + (mi <-> m 3 ) + (mi <-> 7714) 

(B.11) 



All dependence upon the renormalization scale fx cancels out from the total supersymmetric 
matching conditions in y SUSY and Z SUSY . 
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Figure Captions 

Fig. 1. One- loop MSSM (a) bubble, (b) penguin and (c) box graphs which contribute to 
the short distance coefficients of effective theory operators that mediate K L — > 
7r°e+e- and B -> X s £+£- decay. 

Fig. 2. Fractional suppression p(K L — > n°e + e~) of the direct CP violating component 
of the K L — > 7r°e + e _ branching fraction in the MSSM relative to the Stan- 
dard Model prediction. The suppression factor is plotted in (a) as a function of 
Ao(Mqut) and mo(McuT) with m^(m z ) = 90 GeV and tan/? = 5 held fixed. 
The same quantity is displayed in (b) as a function of m^(m z ) and tan/3 with 
Ao(Mgut) = -500 GeV and m (M GVT ) = 100 GeV. Excluded MSSM param- 
eter space points are indicated in both plots by LEGO blocks saturated at their 
maximum values. 

Fig. 3. Ratio of the nonresonant electron asymmetry observable A(B — * X s e + e~)-NR 
calculated in the MSSM relative to its Standard Model value. The ratio is plotted 
as a function of Aq(Mgut) andmo(McuT) with m^(m z ) = 150 GeV and tan f3 = 
20 held fixed. Excluded MSSM parameter space points are indicated in the plot 
by LEGO blocks saturated at their maximum values. 
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